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Abstract 

An interaction of electromagnetic field with metamaterial thin film containing 
split-ring resonators with Josephson junctions is considered. It is shown that 
dynamical self-inductance in a split rings results in reduction of magnetic 
flux through a ring and this reduction is proportional to a time derivative 
of split ring magnetization. Evolution of thin film magnetization taking into 
account dynamical self-inductance is studied. New mechanism for excitation 
of waves in one dimensional array of split-ring resonators with Josephson 
junctions is proposed. Nonlinear magnetic susceptibility of such thin films is 
obtained in the weak amplitude approximation. 

Key words: Metamaterial, Josephson junction, self- inductance, split-ring 
resonator, thin film 
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1. Introduction 

During the last decade metamaterials have become the focus of intensive re- 
search [H, Q, S]- The nonlinear electrodynamics of metamaterials is of special 
interest due to thepresence of new nonlinear phenomena which are specific to 
metamaterials @, H, 0] and due to the fact that, in several cases, nonlinearity 
is a characteristic feature of nanoscale systems 0, S]- The nonlinear response 
of metamaterials can also be the result of deliberate design. For example, 
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nonlinear dielectrics or diodes can be inserted into split rings 0, [h|. Joseph- 
son junctions (JJ) are known to be strongly nonlinear ll[ with low losses and 
therefore they are a natural way to introduce nonlinearity in metamaterials. 
Electromagnetic field interaction with metamaterials containing split-rings 
with Josephson junctions inserted into the gap was recently considered in 
several papers. For example, localized oscillations in chains and two dimen- 
sional arrays of such split-rings were considered in 12|, [l3|, ll4j . The existence 
of metastable states in a medium of this type and transitions between these 
states was discussed in (l5| . 

Currently metamaterials primarily exist in the form of thin films. Thus, it 
is a clear choice to study the interaction of an electromagnetic wave with thin 
films containing strongly nonlinear structural units. Split-rings containing 
Josephson junctions are natural building blocks for such thin films. This 
subject is considered in the first section of this paper. In particular, we 
investigated the effect of dynamical self-inductance in split rings which results 
in reduction of magnetic flux through a ring. This reduction is proportional to 
the time derivative of split ring magnetization. This effect manifests itself as 
non Fresnellian reflection from a film and represents an additional mechanism 
increasing rate of oscillation relaxation in ring-resonators. We show that this 
additional damping leads to dramatic change in the evolution dynamics of 
film magnetization. 

The interaction of a chain of split-rings containing JJ with an electromag- 
netic field is considered in the second section. We expressed the individual 
current of a particular ring in terms of magnetic fluxes through all other rings 
in a chain. In continuous limit, taking into account near neighbor interac- 
tions to leading order, we obtained sine-Gordon type of equation describing 
"continuous" chain dynamics. The effect of dynamical self inductance is 
presented in this model as an additional damping term. External force in 
this equation contains the second spatial derivatives of the incident field and 
therefore suggests a new mechanism for the excitation of longitudinal oscil- 
lations by a normally incident field (without tangential component along the 
chain) . 

In the third section, we obtained the nonlinear magnetic susceptibility of 
a thin film containing split-rings with JJ in the limit of weak amplitudes. 
This susceptibility describes third harmonic generation and self modulation 
due to high frequency Kerr effect. 



2 



2. Electromagnetic wave interaction with thin films containing Joseph- 
son junctions 

Electromagnetic field interaction with thin films (Ax C A) is a well studied 
subject in the literature 0, H H H H Q Ej ■ Most of these works 



consider thin films being polarized under an external field. In this paper we 
study magnetoactive thin films. This requires the derivation of new modeling 
equations which we present in the following subsection. 

2.1. Transmission and reflection of electromagnetic wave on magnetoactive 
thin films: basic equations 
We consider a plane electromagnetic wave normally incident from — oo 
along the x axis on a magnetoactive thin film. The corresponding Maxwell 
equations have following form: 

dHy = ldD ± dE ± = ldB^ 
dx c dt ' dx c dt 
We take into account only magnetic response. Total magnetic inductance 
reads as 

B y (x,t) = B host (x,t)+±irM(t)6(x). 

Here magnetization of the thin film »s IM^ S ' can be represented as a product 
of surface magnetization AfW and the width of the film /. B host (x,t) is the 
magnetic inductance of the host material. Integrating the equations ([I]) over 
x from — 5x to +5x and taking the limit Sx — > leads to the boundary 
conditions at the point x = 0: 

H y (0-) = Hy(0+), (2) 

E z (0-) - E z (0+) = - — (3) 

c at 

In a homogeneous medium for x < and at x > the system of equa- 
tions (jTJ) can be reduced to the wave equation 

d 2 H _ e_cPE_ 

1h? ~ c~cW ^ ' 

where e — E\ for x < and e = £2 for x > 0. Introducing variables qj = 
k 0y/ /Ej, j — 1,2 and k = lu/c we can represent the solution of (jlj) in terms 
of Fourier components as follows: 

T ~ T , \ \ Ae iqix + Be~ iqix , x<0, , r , 
H ^ = \Ce^, x>0, (5) 
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For the electric field we have 

~ f -fcfafco)- 1 {Ae*** - Be-***) , x < 0, 

1 -ftCeaifcb)- 1 ^*, *>0, (6) 

Taking into account boundary conditions (|2J) and ([3]) we obtain relations 

A + B = C, 

**(A-B) = ^C + 4nk 2 lM^\ 
£l £2 

which define the amplitudes B, C via the amplitude of the incident wave 
A = H in : 

jj(tr) = c= 2 ^qi jjjin) + 4WT^£1£2 ^ {a) 

£2qi + £iq2 £2qi + £\q2 

e 2 qi + £iq 2 £ 2 qi + £i<?2 

In the simplest case when Si = e 2 = £ and dispersion of the host material 
is negligible (e = const), expression §J§ can be rewritten as 



C 

In temporal and spatial variables it reads as 

H , {t) = ff(i „> (() _ Mv^aMO (g) 

The expression for magnetic field inside the film H(t) = follows from the 
continuity condition for the tangential components of the magnetic field (j2J). 
For further analysis we need to determine properties of the magnetization 
for the film containing Josephson junctions. 

2.2. Magnetic response of split-rings with Josephson junctions 

We consider a thin film composed of split rings with Josephson junctions 
in their gaps (see Fig. [1]). Orientation of the magnetic field is orthogonal to 
the split ring's plane. Fig. [2] shows an equivalent electric circuit of a split 



ring with Josephson junction [25|]. Here £ stands for electromotive force, L 



is inductance of a split ring, C and R are capacitance and resistance of a 
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H 



Figure 1: Schematic of a split ring with Josepson junction in the gap. Magnetic field H is 
orthogonal to the plane of a split ring. 

L 




Figure 2: Equivalent electric circuit of a ring with Josephson junction. Here E stands for 
electromotive force, L is inductance of a split ring, C and R are capacitance and resistance 
of a Josephson junction. 



Josephson junction. The electromotive force connected with the magnetic 
flux $ is in accordance with Faraday's law: 

£ = -. 10 

c dt v ; 

Self-induced electromotive force inductance Ul can be expressed via current 
J through inductance L: 

Josephson voltage Uc and current Jj are defined by following expressions 
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h dcf) 

Uc= 2e~dl' J, = Jo sin (12) 
In accordance to Kirchhoff laws the sum of the electrical potential differences 
in a circuit is equal to an electromotive force 

£ = Ul + U c , (13) 

and the sum of currents flowing towards that point is equal to the sum of 
currents flowing away from that point, i.e. 

Jj + J R + J C = J- (14) 
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The equation ([131) gives 



|* + ^+s* = ° (15) 



As follows from (1151) the current J can be expressed in the following form 

J=-£ (16) 

Let us consider the second Kirchhoff equation (|14|) . Variables in the equa- 
tion (CHI) are defined as follows 

Jc = c dUc - c hd2(t) 



dt 2e dt 2 ' 
Uc h d(p 
R = ~R = 2e~R~ctt 

The second Kirchhoff equation (TH|) transforms to 

c? 2 ^ 1 d(b c 2 , 2eJ . , 2ec . 
^ + RCTt + LC^ + ~W Sm = " ftLC* (17) 

Magnetic flux is determined by the external magnetic field H and by a cross 
section of a split ring with radius a: 



$ = / Hds « 7ra 2 #. (18) 
is 1 

The modulus of magnetization vector can be defined as 

^a 2 , , 

M = p J, 19 

c 

where p is the density of currents (or contours). Thus magnetic response is 
governed by the following system 

d 2 (j) 1 d<p c 2 2eJ , 2ec7ra 2 rr , 

dP + iVdi + W* + lv™* = -l[LC H (20) 

7ra 2 p f he \ , , 

" = — r (21) 
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Flux and magnetic field in (|2"U|) and (121 j) can be expressed via the external 
magnetic field using equation (JHJ). Thus equations (1201) and (|2"Tj) can be 
rewritten in the following form 



^2 



^ 2 



M (s) 



T^rr + oj^cf) + $ sin i 



~dt 

2ec(na 2 ) 
CLh 



dt 



ira 2 p 



he 

Ye' 



+ ira 



(22) 
(23) 



Here ujt is Thomson frequency (uj, = c 2 /CL), coefficient T = 1/CR describes 
dissipation, and $ = (2eJo)/(Ch) determines strength of nonlinearity. The 
coefficient h/2e can be represented via the quantum of magnetic flux $ = 
nhc/e. Let us introduce dimensionless variables using following scaling: 



H in (t) 
M(t) = 



he 



2eira 2 
hciia 2 p 
2eL 



h(t/r ), 
m(t/r ). 



T 



L), 



T > 



The system of equations fl22|) and (123]) in new variables reads 

d 2( t> d(f) , . , ., . 
+ + <P + K>sm(f) = - [h - dip) , 



i> = - 



here r = u^t, 7 = Y/ujt, k 



d<p dh 
dr dr 



•&Ilo\ and 



(24) 



5 = 1 



Any/ep{iia 2 ) 2 tjj T 



cL 



I 



47rv / ep(vra 2 ) 2 
Ly/LC 



This model describes the magnetic response of a thin film with a diluted 
concentration of split rings containing Josephson junctions. Interaction of 
split rings in this case occurs only via the external electromagnetic field. 
Effects of near neighbor interaction between split rings were considered in 



the literature and can be found in [L2j, [13 
in 



28|, |2<J, |14J. Work presented 



121 . 1281 . 129j considers interacting split rings without Josephson junctions. 



Dense arrangement of split rings with Josephson junctions, which requires 
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consideration of near neighbor interactions is presented in 13|, [14| . In these 
papers the external force acting on an array of split rings was determined only 
by external magnetic field. Based on this work we derive a generalization 



of the model presented in [13|, [14|. Our equations take into account the 
influence of an additional magnetic field due to induced magnetization of a 
split rings. This additional effect is accounted for by the second term in the 
equation Q. This effect results in an increase of the energy dissipation rate 
due to field radiation from the film. An additional dissipation changes the 
magnetization relaxation process to the equilibrium states. Additionally, the 
effect of induced magnetization results in a strong dependance of relaxation 
process on the frequency of external field. 

2. 3. Impact of dynamical magnetic self inductance 

Evolution of a magnetic field described by the equations can be 
illustrated in terms of the dynamics of a Newtonian particle in the potential 



u(4>) = ^ 2 



K COS< 



This potential can have different numbers of minima (stable stationary points). 
The number of such minima is determined by the value of k. Fig. [3] illus- 
trates a potential with one minimum (solid line), where k = 1, and with three 
minima (dashed line), where k = 8. Two maxima corresponding to the last 
case describe unstable states. First we consider potential with one minimum. 
From Equations (1241) it follows that taking into account the effect of self in- 

U 




Figure 3: Graph illustrates the dependance of potential U((f>) — <fi 2 /2 — kcos0 on <fr. The 
potential can have several minima. Solid line illustrates one minimum potential (k = 1), 
dashed line illustrates three minimum potential (k = 8). 

duced magnetization results in the increase of damping factor 7 — > 7+<5. This 
effect is illustrated in Fig.Hl The left subfigure in Fig. H] shows the relaxation 
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of magnetization without taking into account self induced magnetization, and 
the right subfigure shows relaxation in the presence of self induced magneti- 
zation. The relaxation process, which is shown in the right subfigure is faster 
than relaxation in the left figure. Magnetization in this case is excited by an 
external magnetic field of a gaussian shape h(t) = aexp [— (i/io)]> a — 4.5 and 
to = 0.5. Parameters of the equations (I24p for this example have been chosen 
as follows: k — 1, 7 = 0.05, 5 = 0.02 Second, we consider the evolution 




Figure 4: Magnetization relaxation in the case when self induced magnetization of a 
split ring is not taken into account - left figure and is taken into account -right figure. 
Magnetization is excited by by incident spike of the gaussian shape h(t) = aexp [— (t/to)], 
a = 4.5 and £o = 0.5. Evolution takes place in a single minimum potential U{<j)) = 
(j> 2 /2 — cos0 (k = 1 - solid line on Fig. [3]). 

of magnetization in the case of a three minimum potential. The difference 
in the evolution of magnetization between the two cases with and without 
self induced magnetization is more dramatic in a three minimum potential. 
There are four oscillatory regimes corresponding to oscillations around the 
central minimum, around the two side minima and when oscillations are tak- 
ing place above two local maxima of the potential (see Fig. [3]). These types 
of oscillations have different frequencies and, in the last case, the period of 
oscillations is largest and the amplitude is limited from below. Fig. shows 
the dynamics of magnetization with and without self induced magnetization 
for the same values of parameters and external magnetic field. The initially 
excited oscillations corresponding to the trajectories of Newtonian particles 
above local minima of the potential are decaying and eventually switching 
to faster and lower amplitude oscillations around one of the minima. The 
moment of switching and asymptotic stable state is highly sensitive to the 
presence of self induced magnetization. In this particular case, which is illus- 
trated in Fig. [5], switching to oscillations around the central minimum occurs 
after seven cycles of large amplitude oscillations - left subfigure. The pres- 
ence of an additional decay in relaxation of magnetization due to self induced 
magnetization, which results in excessive radiation, results in the switching 
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Figure 5: Magnetization relaxation in the case when self induced magnetization of a 
split ring is not taken into account - left figure and is taken into account -right figure. 
Magnetization is excited by by incident spike of the gaussian shape hit) = aexp [— (t/to)], 
a = 12.1 and to = 0.5. Evolution takes place in a two minima potential U(<f) = 4> 2 /2 — 
k cos (f) (k — 8 - dashed line on Fig. [3]). 



to oscillations around right minima only after three and a half oscillatory 
cycles - right subfigure. In both cases oscillations are excited by an external 
spike of the magnetic field h(t) = aexp [—(t/to)], a, = 12.1 and to = 0.5, 
where k = 8. 



2.4- Multistability of magnetization 



The mirrorless optical bistability has been predicted in a bulk [18| and in 



the thin film [19( of dipole-dipole interacting two-level atoms. We can expect 
this effect in the case under consideration too. In the case of stationary fields, 
system ( 12^1) reads as 

(m + h) + /tsin(m + h) = h, h = const, (25) 

This is an implicit equation for <p as function of h. The general solution 
of ff25l) can have several roots. For a three minimum potential, the solution is 
illustrated in Fig.O The presence of three stable and two unstable stationary 
solutions leads to the multivalued nature of the function 4>(h), which results in 
multistability of magnetization. This function consists of stable and unstable 
branches, which lead to hierarchical set of hysteresis loops. The simplest first 
order hysteresis loop is shown in Fig. [6j 



3. Near neighbor interaction interaction between split rings 



Following the work in 12] let us take into account the magnetic dipole- 
dipole interaction between neighboring SRRs characterized by cross-inductance 
Lcross- Note that cross- inductance decays as the cube of the distance. Let us 
restrict ourself to only nearest neighbor SRR interactions. 
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10 20 



Figure 6: Graph illustrates multistability of magnetization m versus value of stationary 
magnetic field h. Parameters of three minima potential arc k — 8 - dashed line on Fig. [3]). 



The voltage on the nth junction is expressed via phase 4> n as 

f/ nc = A| (26) 
' 2e dt K ' 

The voltage on the junction is equal to U n> c and is also equal to the voltage 
U Ut R. Thus the normal (i.e., non superfluid) current thought the junction is 

U n , c /R- 

According to Kirchhoff law, voltage on the ring is equal to the sum of 
acting electromotive forces. Therefore, the first Kirchhoff's equation has 
following form: 

^ = _ l ^ _ l_ ^ + aj^y Zn (27) 

Here £ n = —c~ x d§ n jdt is electromotive force induction generated by high 
frequency magnetic flux acting on nth ring, I n ±\ is the current in the (n±l)th 
ring. 

The second Kirchhoff equation corresponds to charge conservation: 

In = Jn + C n ^ + ^ (28) 

at R 

Combining equation ([2611 with (I2TI) we obtain 

c ( cH \ 

In + V (In-1 + In+l) =~J j [^n + ~ <\>n J , (29) 
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where r\ = L cross /L. This equation must be complemented with equations 
defining currents at each split-ring: 



Chd 2 (p n h d(j) r 
2e dt 2 + 2er dt 



I n = Jon sin (j) n + - — ^ + -P-. (30) 



Expression fl29|) represents a linear recurrency relation. In the limit TV 
oo, we can employ the generating function method. 

3.1. Calculation of individual currents in a split rings 
Let us introduce notation 

c / _ ch 



L \ 2e 
and rewrite equation ( 1291) as 

I n + 7] (7 n _! + 7 n+ i) = F n , (31) 
Introducing the generating function 

+00 

P(y) = V"' 

n=— oo 

defining the function 

+oo 

F(y) = F ne iny (32) 

n=— oo 

and using equation f[3"Tj) we can find relation between P(y) and 

P( y ) = EM . (33) 

Kyj l + 2r]cosy K ' 

This expression can be represented as a power series 

oo 
m=0 
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Function cos m (y) can be expressed as 

cos m {y) = 2 - m e imy (1 + e- 2iy ) m = 

= 2~ m e imV l m \ e ~2ipy 

then the generating function has the form: 



co m s \ 

p( y ) = j2(-^rJ2{ m ) e " lv{m ~ 2p)F ^- 

m=0 p=0 \P ' 



P= 

The current in the n-th ring is determined by the formula 

/+oo 
P(y)e my dy. 
-co 

Using the expression for F(y) fl32|) we can find that 



co m r \ 

m=0 p=0 \ ^ ' 



(*)■ (34) 



The expression for I n can be simplified taking into account smallness of the 
parameter 77 <§C 1. Up to the second order of 77, the expression for I n (t) has 
the form 

I n (t) « F n (t)-r ? (F n+1 (t)+F n _ 1 (t)) + 

+ r? 2 (F n+2 (t) + 2F ft (t) + F„_ 2 (t)) . (35) 

5.^. Continuum chain approximation 

In the continuum approximation we can introduce a new variable y = nl r , 
where l r is the distance between neighboring rings, and represent F n ±\(t) as 

dF l 2 d 2 F 

Limiting ourself to the first order expansion with respect to rj we obtain from 
(J35D 

I(y,t) « (1- 2 V )F(y, t) - V l 2 -^F(y,t) 



dy 2 

F(y,t)-r ] l 2 r ^-F(y,t). (36) 



' dy 2 ' 
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In our model F(y,t) was denned as 

F(y,t) = ~ U(y,t) + ^<j>(y,t) 

Since the magnetic flux through ring is greater than the flux through the 
junction (gap), we can use the approximate expression $(?/, t) « (na 2 )H(y, t). 
Substitution of this expression in f[3"Uj) leads to 

j . , Chd 2 4> n h d<j> n 
J 0n sin <p n + — — + 



2e <9t 2 2eR dt 

~[Hy,t) + ^4>(y,t)) + 



L \ dy 2 2e dy 

Using dimensionless variables defined above, the resulting system of equa- 
tions can be written as 

d 2 4> d 2 4> d(j) . . 

+ 7t; — h (ft + K sm (ft = 



dr 2 <9£ 2 Or 

= -(h-5^) + — (h-5ij), (38) 
fd<j> dh\ , , 

Here we introduce dimensionless spatial variable £ = y/(l r y/v)- 

In case of plane electromagnetic wave propagating along x, dependence 
of h on y vanishes. Spatial derivatives in the equation (1391) disappear and 
equation describes homogeneous dynamics of the chain as function of time. 
Temporal-spatial chain dynamics takes place when the chain is excited by an 
electromagnetic beam with spatial distribution of intensity across the beam. 
In the general case, excitation of surface waves requires matching of the 
surface wave vector with tangential component of the vector of the incident 
wave. Equation ( 1381) demonstrates that sharp spatial gradients across the 
beam act as an external force for oscillations along the chain. 

4. Small-amplitude approximation in SRRs 

In the linear approximation equation (1221) reads 



d 2 (h 1 dd> ~ , 2ec ^ . . . . 



14 



where Q is renormalized Thompson frequency Q 2 = (LC/c 2 ) 1 + (2eJo/Ch). 
The Fourier components of current are represented as 



na 2 c ( ujn 



Fourier components of magnetic induction B = H + 4irM are expressed 
through magnetic permittivity B = [i(cu)H, which leads to 



fi(u) = 1 + 4im ri 



(na 2 ) 2 I ^-1 



L \Q? — uj 2 — iToo 



Since equation (1221) is nonlinear, the magnetic response is described by a 
nonlinear coupling between field H and magnetization M. The total magne- 
tization can be represented as sum of a linear part Mn n and nonlinear part 
M n \. Taking into account the definition of magnetization (1231) we obtain 
following expressions 



Tia 2 ( o TT he 
M hn = -n m — \na H + — 

M nl = -n m — — 
L 2e 

where ft ^ is solution of the linear equation, and ft 1 ' is correction to solution 
of the linear equation (|40l) . 

Expansion of sine-function up to cubic term 3 ), transforms (j22j) into the 
Duffing equation: 

% + + + * -ZD?® («, 

where = 2eJo/Qch is parameter characterizing nonlinear response of Duffing 
oscillator. Let us consider only the first correction term for solution of (122]) . 
Substituting <\> = ft ^ + ft 1 ^ into (|22p and collecting terms of same order we 
obtain: 

Qifto) _ ^(o) _ 2ec 



a* +r ir + ^ --ozn^ (42) 



at 2 at 



+ r ^ + n 2 ft^ = -#^ 3 . (43) 
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Solution of (j42p gives a zero order approximation for (ft: 

( />« s )(u) = -5 m £(u)H(u), (44) 

where 

2ircea 2 liuJL . 2 \ 
is a coupling parameter. Function 

1 



f2 2 — u) 2 — icoT 



is the standard Lorentzian function. For harmonic wave with carrier fre- 
quency ujq, we have 

H(uj) = H 8(u - ujq) + Hq5{uj + u ). 

In this case 

0(0) (a,) = _«j m [£(uj)H 5(uj - w ) + £(-w)/7*5(w + u )\ , 

and 



(t) = -6 m [^(uo)H e-^ ot + £(-uj )H*e iuJot ] . (45) 

The next step is to substitute expression for <f>(°'{t) into the right hand 
part of the equation ( 1431) : 

^0( ) 3 (t) = a^Hfe- 3 ™ + 

+ a 2 (u; )\H \ 2 H e- iuJ0 + c.c. 

Here oti(wo) and 62(^0) are defined as follows: 

a 2 (uj ) = 3^£(w )£(w )£(-w ). 

From (J4"3l) . taking account of the results obtained above, we can find 
0«(w) 

0«(w) = -a 1 (w )£(w )if 3 5(w-3wo)- 

-a 2 (uJo)£(u}o)\H \ 2 H 5(uj - u ) + c.c. 
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Having (j)(°\oj) and <j)^(co) we can write expressions for the linear and 
nonlinear parts of magnetization: 




Here we use notations for nonlinear magnetic susceptibilities of third order 



The first term in the expression for M n i(uo) describes process of third har- 
monic generation, second term in this expression describes phenomena of 
Kerr's self-modulation. 

5. Conclusion 

We derived equations describing electromagnetic pulse interaction with 
thin films containing split rings with Joshephson junctions, taking into ac- 
count dynamical self-inductance in the split rings, and analyzed the impact 
of dynamical self-inductance on the transmitted and reflected wave. We 
demonstrated an increase of magnetization relaxation rate due to dynamical 
self-inductance. This additional damping can results in significant change 
in the evolution dynamics of film magnetization. We also derived an equa- 
tion describing the interaction of a chain of split-rings containing JJ with 
an electromagnetic field and showed that the current in a particular ring is 
defined by the magnetic fluxes through all other rings in a chain. The con- 
tinuous limit transforms this model into sin-Gordon type of equation. This 
equation is driven by an external force that contains second spatial deriva- 
tives of the incident field. Presence of the field spatial derivatives suggests 
a mechanism for the excitation of magnetization waves along the chain by 
normally incident field. We found analytic expression for nonlinear magnetic 
susceptibility of such films in the limit of weak amplitudes, describing third 
harmonic generation and self modulation due to high frequency Kerr effect. 



a 2 = oi 




W(3uj;u ,u; ,u; ) = #5 3 m [£(u )f £(3^ ), 
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